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Abstract: 

The pseudo SU(3) approach is used to describe low lying states and BE(2) intensities of 
rare earth and actinide nuclei which are (3(3 decay candidates. The (3(3 half lives of some of 
these nuclei to the ground and excited states of the final ones are evaluated for both the two 
and zero neutrino emitting modes. The existence of selection rules which strongly restricts the 
decays is discussed. These restrictions represent a possible test of the model. Up to now the 
predictions are in good agreement with the available experimental data. 

Resumen: 

El formalismo pseudo SU(3) es usado para describir los estados de menor energia y sus 
intensidades BE(2) de micleos en la region de las tierras raras y los actinfdos que son candidatos 
a decaimientos (3(3. Las vidas medias de esta desintegracion para algunos de estos micleos al 
estado fundamental y estados excitados del nucleo hijo son evaluadas en los modos con y sin 
emision de neutrinos. Se discute la existencia de reglas de seleccion que restringen fuertemente 
los decaimientos y representan una posible prueba para el modelo. Hasta el momento el acuerdo 
entre las predicciones y los datos experimentales disponibles es bueno. 

PACS: 23.40.Hc, 21.60.Fw 
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1 Introduction 



The detection of the neutrinoless double beta decay (PPou) would imply an indisputable evidence 
of physics beyond the standard model and would be useful in order to select Grand Unification 
Theories [|l]]. Theoretical nuclear matrix elements are needed to convert experimental half-life 
limits, which are available for many /?/5-unstable isotopes0, into constrains for particle physics 
parameters such as the effective Majorana mass of the neutrino and the contribution of right- 
handed currents to the weak interaction. The two neutrino mode of the double beta decay 
{PP2u) is allowed as a second order process in the standard model. It has been detected in nine 
nuclei and has served as a test of a variety of nuclear models. It is the best available proof 
we can impose to a nuclear model used to predict the PPou matrix elements. 

Many experimental groups have reported measurements of (3l32u processes[0]. Nearly for all 
the cases the ground state (g.s.) to ground state (0+ 0^) decay was investigated. In direct- 
counting experiments the analysis of the sum-energy spectrum of the emitted electrons allowed 
the identification of the different /3/3-decay modes. The experimental research on the double 
beta decay of heavy deformed nuclei enjoys an increasing activity. Thus, the l3l32u half-life of 
and ^^''A^(i[^ ^ Q have been measured and search is in progress for ^'*^P'u|^. 

In previous papers 0, H we used the pseudo SU(3) shell model to evaluate the two neutrino 
double beta half lifes of eleven heavy deformed potential double beta emitters. We found good 
agreement with the available experimental information. The radiochemically measured ^38^ 
decay raised some expectations, given the experimental arrangement cannot discriminate 
between the different double beta decay modes, and some theoretical estimates^] predicted 
similar decay ratios for both the Oz/ and the 2z/ modes. Our calculations for the two neutrino 
mode were consistent with the experimental result. 

We have also calculated the zero neutrino matrix elements for six heavy deformed double 



beta emitters||T^. Using the upper limit for the neutrino mass we also estimated their (3(3ou half- 
lives. In the case of ^ssf/ found the zero neutrino half-life at least three orders of magnitude 
greater than the two neutrino one, giving strong support of the identification of the observed 
half-life as being the two neutrino double beta decay. The (3/32^ of ^^^Nd to the ground and 



2 



excited states of ^^^Sm has also been studied|TT[|. In the present paper we review the above 
mentioned calculations as well as the pseudo SU(3) formalism used in those works. 

2 The pseudo SU(3) formalism 

In the pseudo SU (3) shell model coupling scheme[0] normal parity orbitals [r], I, j) are identified 
with orbitals of a harmonic oscillator of one quanta less fj = t] — 1. This set of orbitals with 
j = j = I + s, pseudo spin 5 = 1/2 and pseudo orbital angular momentum / define the so 
called pseudo space. The orbitals with j = / ± 1/2 are nearly degenerate. For configurations of 
identical particles occupying a single j orbital of abnormal parity a convenient characterization 
of states is made by means of the seniority coupling scheme. 

The many particle states of Ua nucleons in a given shell rja, a = z/ or tt, can be defined by 
the totally antisymmetric irreducible representations {1"" } and {1"°} of unitary groups. The 
dimensions of the normal (A^) parity space is = {fja + l)(??a + 2) and that of the unique {A) 
space is fi^ = 2170 + 4 with the constraint = n^ + n^. Proton and neutron states are coupled 
to angular momentum and in both the normal and unique parity sectors, respectively. 
The wave function of the many-particle state with angular momentum J and projection M is 
expressed as a direct product of the normal and unique parity ones, as: 

\jM>= y: >(E)\j^ >Vm (1) 

We are interested in the low energy states which have J = 0, 2, 3, 4, 6. In the normal subspace 
only pseudo spin zero configurations are taken into account. Additionally in the abnormal parity 
space only seniority zero configurations are taken into account. This simplification implies that 
J:^ = = 0. This is a very strong assumption quite useful in order to simplify the calculations. 
Its effects upon the present calculation are discussed below. 

The double beta decay, when described in the pseudo SU(3) scheme, is strongly depen- 
dent on the occupation numbers for protons and neutrons in the normal and abnormal parity 
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states ,n^,n^^. These numbers are determined filling the Nilsson levels from below, 

as discussed in p. In particular the PP2U decay is allowed only if they fulfil the following 
relationships 

</ = <^ + 2 , = . (n^ 

If they do not, the decay becomes forbidden. This is the first selection rule the pseudo 
SU(3) formalism impose to the double beta decay. As an example, for ^^^Nd we have obtained 
the occupation numbers = 4, = 6, = 2, = 6. 

We have selected the standard version of the pseudo SU(3) Hamiltonian P3| . It is constructed 



by a spherical Nilsson hamiltonian which describes the single-particle motion of neutrons or 
protons, a quadrupole-quadrupole interaction and a residual force. The latter allows the fine 
tuning of low lying spectral features like K band splitting and the effective moments of inertia. 

With the occupation numbers and the hamiltonian discussed above, the wave function of 
the deformed ground state of -^^^Nd can be written^] 

|i50^^^ 0+) = I ihu2)t, = = 0; {zrs/2)l, = M^^ = > 



(3) 



{16}^{23}^(12,0)^; {16}423},(18,0),; 1(30, 0)ir = IJ = M = > 



3 The double beta decay 



The inverse half life of the two neutrino mode of the /3/3-decay can be expressed in the form 14 



= GM) I MM) |2 

where G2v{J^) are kinematical factors. They depend on Ej.,j = ^[ 
is the half of total energy released. The nuclear matrix element is evaluated using the pseudo 
SU(3) formalism. For the ^^^Nd^ ^^^Sm case it can be written as 



(4) 

E{Ja)] + rrieC^ which 
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Mf/(J+) = a{J) h{n^) £jy^'^ 

E < mni, {Ofi)ll\\{Xofio)KoJ >i E < (30,0)1 0, {Xofio)KoJ\\{Xfi)alJ >p 

(AoMo)i^o P 



E 

p' 



(12,0) (0,0) (12,0) 1 

(18,0) (Ao/io) (12,2) p' 

(30,0) (Ao/Uo) (AAi)a p 

1 1 1 



<(12,2)|||[ao,)iao,w](^°^»Mll(18,0) 



(5) 



In the above formula < ..,..||,,> denotes the SU(3) Clebsch-Gordan coefficients [p!5|, the 
symbol [...] represents a 9 — A/i recoupling coefficient , < .. 1 1 1 .. 1 1 1 .. > is the triple reduced 
matrix elements [jl^ and 



a(0) 



4^ 

(2r)+l),/2^FT' 



a(2) 



2r,+l 



5r,(7?-l))2r?-3) 
3(2»y+l) 



1/2 



K^;:') = [« + 2)(r/ + i-n;;^/2)]V2, 



c 



0.70 

A173 



[2Z+ 1-0.76((Z + 1)4/3 _^4/3)]j^^g|^ _ 



(6) 

The SU(3) tensorial componentes (Ao,/io) of fhs normal part of the double Gamow- Teller op- 
erator must be able to couple the proton and total irreps (18,0) and (30,0) associated with the 
ground state of ^^^Nd to the corresponding irreps (12,2) and (A/i)(j) which characterize the 
ground and excited rotational bands in ^^^Sm. If it is not possible the (3(32^ decay to the O"*" 
and 2+ states of an excited band is forbidden. This is the second selection rule imposed by the 
model to the decay. 

In the case the Oz/ decay exists, the virtual neutrino must be emitted in one vertex, and 
absorbed in the other. Since in the standard theory the emitted particle is a right-handed 
antineutrino and the absorbed one a left-handed neutrino the process requires that a) the 
exchanged neutrino is a Majorana particle and b) both neutrinos have a common helicity 
component. The helicity matching can be satisfied in two ways: a) the neutrinos have a 
nonvanishing mass and therefore a "wrong" helicity component proportional to m^/E^. The 
decay rate will be proportional to < >^. Or b) the helicity restriction could be satisfied if 
there is a right handed current interaction. In this case a nonvanishing mass allowing mixing 
of neutrino types is also required [0, . 
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The (3(3oi, matrix elements were evaluated using the pseudo SU(3) formalism ||1 1||. The final 
expressions are similar to those related with the two neutrino mode, given in Eq. (5), but 
restricted to the decay to the ground state of the final nuclei. They explicitly include the 
"neutrino potential" originated by the interchange of a Majorana neutrino between the two 
emitting neutrons. The Oz/ mode allows transitions with Al ^ 0, whose contribution to the 
/3/?0!^ matrix elements were found an order of magnitude lesser than those which do not change 
the orbital angular momentum. 

The five nuclei which were found to have the (5(52v mode strictly forbidden within the pseudo 
SU(3) model are allowed to have /?/3o!/ decays because the neutrino potential allows parity 
mixing transitions. We would expect that these Oz/ matrix elements, related with A/ = 1, 3, 5 
transitions were of the same order of magnitude of the A/ = 2,4,6 transitions and about a 
factor ten lesser than the A/ = matrix elements. This hindrance is comparable with the 
modifications that the symplectic model would introduce, allowing more than one shell active 
for each kind of nucleons. These questions are under current investigation. 

4 Results and discussion 



We predict, according to the assumptions of the model, that the two neutrino double beta 
transitions of the following five nuclei: 

















^ 232^ 


244PM - 





are forbidden. In particular we predict a null result for the present search in ^^'^P-u^, a potential 
strong test of our model. Inclusion of pairing within the same truncation scheme will allow 



mixing of different pseudo SU(3) irreps, and is in progress[19|. The lack of a spin-isospin channel 
in our Hamiltonian is perhaps an additional limitation, but it is included, at least partially, in 
the quadrupole-quadrupole interaction, when it is recoupled. 
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In Table 1 are indicated the results for the six jSjS emitters in which we get a j3l32u matrix 
element different from zero^]. The value for the integrated kinematical factors were obtained 
following the procedure indicated by Doi et al|]20| with QA/gv = 1-0. In the second and third 
columns the theoretical and experimental /3/32i/-half lives are given. The agreement with the 
available data for ^^'^Nd and 232f/ jg good. Improvements in the experimental information about 
the decay of these nuclei together with ^'^^Nd is possible, but the other three are almost excluded 
of detection because their lepton kinematical factors G2V are very small. This is particularly 
true for ^^^Nd, which has the smallest Q value for the (3 (3 decay and thus becomes strongly 
inhibited. 

We also evaluated the jSPoi, matrix elements between the the same six nuclei [|l3]. In the last 
two columns of Table 1 the theoretical predictions and experimental lower limits of the jSjSo^- 
half lives are given. In order to calculate the zero- neutrino half-life we assumed < rrii, >= leV, 
which is the larger neutrino mass parameter compatible with the ^^Ge experiment [^. In the 
case of 238f/ ^]^g predicted Oz/ half life is three orders of magnitude greater than the predicted 
2i> half life, which essentially agrees with the experimental one, confirming that the observed 
jSP decay of ^ssf/ j-^g^g ^j-^g neutrino mode. The transition is strongly dominated by 

a pair of neutrons in the normal parity orbital iii/2 decaying into two protons in the unique 
orbital ^13/2; 

In the case of ^^^Nd, the pseudo SU(3) Ou matrix element reported here is about a factor 
four lesser than the QRPA estimations]^. This is a very relevant result. First, it exhibits 
the stability of the neutrinoless double beta decay matrix elements evaluated in quite different 
nuclear models, in the case of deformed nuclei. Second, this factor of four, which is little 
compared with the order of magnitude variations in the 2i> theoretical estimations, is still 
important in order to extract the parameter < rrie >. In this case the transition is dominated 
by a pair of neutrons in the normal parity orbital hg^2 decaying into two protons in the unique 
orbital h'li/2j again resembling the two neutrino case. 

1 /2 

As can be seen in the last two columns of Table 1, the Tq[, predicted for < >= leV are 
at least three order of magnitude greater than the experimental limits. These results reflect 
the fact that at, the present stage of the experimental (3f3 research, the limits < rrii, >< l.leV 
obtained by the Heidelberg-Moscow collaboration ||2T[ using significative volumes of ultrapure 
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''^Ge are the most sensitive. But, if the PPqi, decay is observed in ^^Ge, at least a second 
observation will be essential, and ^^^Nd is a likely candidate to do this job [^, 0. In the next 
few years the limit for < m^, > extracted from PPq^, experiments is expected to be improved up 



to 0.1 eV and ^^^Nd is one of the selected isotopes 



Finally we review the two neutrino mode of the double beta decay l3l32u of ^^^Nd into the 
ground state, the first excited 2+ and the first and second excited 0"*" states of Sm^^^ Sm^l\. 

In Table 2 the matrix elements and predicted half-lives for the l3l32u decay of ^^^Nd to the 
ground state, the first 2+ and the first and second excited O"*" states of ^^^Sm are presented. The 
matrix elements are given in units of (mgC^)"*-"^'^^''. It must be mentioned that the phase space 
factors differ in about 10% with those presented in Table lp| where a different renormalization 
procedure was used. 

The /3/?2i/ decay to the first excited 0"*" state is forbidden. In this model this is imposed by the 
fulfilment of an exact selection rule. The pair of annihilation operators o-^qJ)!, when expanded 
in their SU(3) components, cannot couple the ^^^Nd g.s. irrep (30,0) to the irrep (20,4) which 
we associated with the first excited O"*" state. Thus the transition between members of these 
particular irreps are forbidden. 

The decay to the second excited 0^ state is allowed but strongly cancelled. The predicted 
half-life is four orders of magnitude larger than that of the decay to the g.s. The f3(32u decay 
to the 2"*" state is inhibited by the fi% dependence of the matrix element as it is discussed in 
Section 3. The matrix element of the /?/?2j/decay to the first excited 2+ state M2i/(2+^ ) is three 
orders of magnitude lesser than the matrix element of the decay to the g.s. M2u{0'^g ). 

The present results contradicts those previously published p3|, ^ were it was speculated that 
the /?/?2v decay of ^^^Nd to the first excited 0"*" state of ^^^Sm could have a similar intensity of 
that to the g.s. We found that in the present formalism this decay is forbidden. If we select 
different occupation numbers for both ^^^Nd and ^^^Sm, taken the deformation of the latter 
nucleus instead of that of the former, we found very similar results for the decay to the g.s. 
and the first 2+ state, but the matrix elements of the decay to the first and second 0"*" states 
becomes interchanged with essentially the same values. Considering the difference in the phase 
space integrals we predict a half live of the order 10^^ years for the decay to the first excited 
0^ state and the decay to the second excited one becomes forbidden. 
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The above discussed reduction of the matrix element of the (3l32u decay to the excited O"*" 
state as compared with the decay to the g.s. is not a general result of the pseudo SU(3) 
scheme. A recent analysis of the case of ^"^Moj^^ shows that both matrix elements are very 
similar and that they are in agreement with the experimental information. In conclusion, 
the appearance of selection rules which can produce the suppression of the matrix elements 
governing a I3(32v transition is a consequence of the details of the irreps involved. 

5 Conclusions 

The pseudo SU(3) model is a very powerful machinery to describe the collective behavior of 
heavy deformed nuclei. It has been used to reproduce very accurately the rotational spectra 
of heavy deformed nuclei, including the K-band splitting, the amplitudes for transitions of the 
E2, Ml and M3 type. 

We have used this model to evaluate the [3[32u half lives of 11 heavy deformed nuclei, using 
the pseudo SU(3) approach together with a summation method, which avoids the closure ap- 
proximation, obtaining good agreement with the available data and making testable predictions, 
quite different from the QRPA ones in some cases. 

We have also evaluated the [3[3q„ half lives of six of these heavy deformed nuclei. We obtained 
an order of magnitude agreement with the QRPA estimations, and a factor four of difference. 
We exhibited predictions for the neutrinoless double beta half-lives assuming < niy >= leV, 
and discussed the relevance of our results. In the case of '^^^U this result complements those 
obtained for the two neutrino jS/S decay, and it confirms the observed half-life as two neutrino 
in origin. 

At last we have studied the decay mode of ^^'^Nd to the ground and excited states of 
^^^Sm. The (3(32u decay to the first excited 0"*" state was found forbidden in the model and the 
decay to the second excited O"*" state has a half-life four orders of magnitude greater than that 
to the g.s.. The decay to the 2+ state is strongly inhibited due to the energy dependence of the 
matrix elements M2u{'^^), two powers greater than that of the matrix element M2u{0~^). 

In all the calculations only one active shell was allowed for protons, and one for neutrons. 
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This is a very strong truncation, of the same type as used in shell model calculations. An 
important consequence of this truncation is the fact that only one uncorrelated Gamow- Teller 
transition is allowed: that which removes a neutron from a normal parity state with maximum 
angular momentum, and creates a proton in the intruder shell ^11/2 ^^^^ earth 

nuclei, iiii2 ^13/2 actinides). This unique Gamow- Teller transition controls the double 
beta decay. If the occupation of the Nilsson levels is such that the number of protons in the 
abnormal states does not change going from the initial to the final state configurations, the 
decay becomes forbidden. 

The pseudo SU(3) model uses a quite restrictive Hilbert space. The model could be improved 
by incorporating mixing between different irreps, via pairing by example ||19||. Also other active 
shells can be taken into account in the symplectic extension ||2^. In both cases the selection 
rules that impose such strong restrictions to the decays of some nuclei can be superseded. 
However if the main part of the wave function is well represented by the pseudo SU(3) model 
those forbidden decays will have, in the better case, matrix elements that will be no greater 
than 20% of the allowed ones, resulting in at least one order of magnitude cancellation in the 
half-life. 
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Table Captions 

Table 1. Theoretical estimates for the /3/3-decay half-hfe in the 2v and mode for several 
heavy deformed nuclei are given and compared with the available experimental data. 

Table 2 The dimensionless matrix elements and predicted half-lives for the /3/32i/ decay of 
^^^Nd to the ground state, the first 2+ and the first and second excited 0+ states of ^^^Sm. 
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Table 1 



Transition 


1/2 

theo 


[yr\ 

exp 


1/2 

theo 


r 1 
[yr\ 

exp 






2.1 X 10^1 




1.18 X 10^8 




''^Nd - 


-^1^8 Sm 


6.0 X 10^0 




6.75 X 10^4 




isoivrf - 


-^150 Sm 


6.0 X lO^s 


9 - 17 X 10^8 i i i 


1.05 X 10^4 


> 2.1 X 10^1 § 




_^186 03 


6.1 X 10^4 




5.13 X 10^5 


> 2.3 X 10^0 26 


1920s - 




9.0 X 10^5 




3.28 X 10^6 




238 _ 


^238 


1.4 X 10^1 


2 X 10^1 [| 


1.03 X 10^4 


> 2.0 X 10^1 II 



Table 2 









0+ - 


-^0^ 


-{g.s.) 


.0549 


6.73 X 10^8 


0+ - 


-^0^ 


-(1) 





OO 


0+ - 


-^0^ 


-(2) 


.00499 


4.31 X 10^2 


0+ - 


-^2^ 




5.38 X 10-5 


7.21 X 10^^ 
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